Random matrices with external source and multiple orthogonal 

polynomials 



(N 



P.M. Bleher and A.B.J. Kuijlaars 
July 25, 2003 

rn ' 
O , 

O , Abstract 

We show that the average characteristic polynomial P„(z) = E[det(z/— A/)] of the random 
Hermitian matrix ensemble exp(— Tr(y(A/) — AMy)dM is characterized by multiple 
orthogonality conditions that depend on the eigenvalues of the external source A. For each 
OO ' eigenvalue aj of there is a weight and P„ has rij orthogonality conditions with respect 

to this weight, if Uj is the multiplicity of aj. The eigenvalue correlation functions have 
determinantal form, as shown by Zinn- Justin. Here we give a different expression for the 
^ . kernel. We derive a Christoffel-Darboux formula in case A has two distinct eigenvalues, 

l/^ ' which leads to a compact formula in terms of a Riemann-Hilbert problem that is satisfied 

. by multiple orthogonal polynomials. 

^ ■ 1 Random matrices with external source 
O 

' Following Brezin and Hikami El [3 El and P. Zinn- Justin |2Ul I21j we consider a random matrix 

p) ^1 ensemble with an external source, 

S '■ J_e-^<V(M)-AM)^M (1.1) 

defined on n x n Hermitian matrices M. The ensemble ()1.1() consists of a general unitary 
invariant part V{M) and an extra term AM where ^ is a fixed n x n Hermitian matrix, the 
^ . external source or the external field. Due to the external source, the ensemble (|1.1() is not unitary 

^ I invariant. For the special Gaussian case V{x) = ^x^, we can write M in (|1.1() as M = H + A 

where H is a random matrix from the GUE ensemble, and A is deterministic, hence in this case 
it reduces to the class of deterministic plus random matrices studied in |^1 ^3 [^1 101 13 IHl IH] ■ 

Zinn-Justin [20] showed that the eigenvalue correlations of ensemble Hl.l|) can be expressed 
in the determinantal form, 



i?m(Ai, . . . , Xm) = det{KniXi, Xj)) 



for some kernel Kn. In this paper, we give a different expression for Kn- We believe that our 
formulation is useful for asymptotic analysis. Indeed, for the Gaussian case V{x) = and 
for the case where A has only two distinct eigenvalues, we have been able to carry out the 
asymptotic analysis almost completely. This will be reported elsewhere. 

Our approach is based on the observation that the average characteristic polynomial 

Pn{z) =E[det(2-M)] 

of the ensemble (|1.1|) can be characterized by the property that 
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for every eigenvalue Uj of A and for k = 0,. . . ,nj — 1, where rij is the multiphcity of aj, see 
Section 2. We can embed the polynomial Pn in a sequence of polynomials {Pk}o where Pk has 
degree k. Then our kernel has the form 

n-l 

Kn{x, y) = e-^(^(-)+^(j')) Pk{x)Qk{y) (1.2) 

fc=0 

where the Qk are certain dual functions (not polynomials in general), see Section 3. 

When A = (no external source), the polynomials Pk are usual monic orthogonal polyno- 
mials with respect to the weig ht e-^(^) on R. In that case, the function is a multiple of Pk 
and the kernel (|1.2|) reduces to the orthogonal polynomial kernel which is familiar in the theory 
of random matrices. By the ChristofFel-Darboux formula we then have 

KJx,y) = c-~¥^M+y(y)) ^'^'^ Pn{x)Pn~l{y) - Pn-l{x)Pn{y) ^ 

In x-y 

where 7^ is the leading coefficient of the orthonormal polynomial of degree k. 

In Section 4 we present an analog of the Christoffel-Darboux formula for the kernel ()1.2() in 
the case where A has only two eigenvalues. We also relate it to a Riemann-Hilbert problem in 
Section 5. 

2 The average characteristic polynomial 

We define the monic polynomial 

Pn{z) =E[det(z-M)] 
where the expectation is with respect to the ensemble 

Proposition 2.1. Suppose A has eigenvalues aj, j = I, . . . ,n with ai ^ aj if i ^ j- Then the 
following hold. 

(a) There is a constant Zn such that 

„ n n 

Pn{z) = — / - ^i) n e-(^("^)-'^^"^)A(A)dA, (2.1) 



where 



A(A) = n(A^ - A, 



and dX = d\id\2 ■ ■ ■ dXn 
(b) Let 

Then we have the determinantal formula 

mio mil ■■■ min 



/oo 
-00 



Pniz) = ^ 



1 



mnO rUnl ■ ■ ■ rUnn 

1 z ••• 



(2.2) 



2 



(c) For j = 1, . . . ,n, 

/oo 
Pn(x)e-(^(^)-'^^-^)(ix = 0, 
-oo 

and these equations uniquely determine the monic polynomial Pn- 



(2.3) 



Proof. Write M = UAU* where U is unitary and A is diagonal, A = diag(Ai, . . . , A„). of M 
on the diagonal. Then by the Weyl integration formula, see e.g. I17j . we have for every 
integrable function / on the space of Hermitian n x n matrices. 



J f{M)dM = 7r-"("-i)/2 l^-Q J I' f(^uAU*)A{XfdXdU 
where dU denotes the normalized Haar measure on the unitary group U{n). Thus 



(2.4) 



Pn{z) 



\{{z - A,) n e'''^'^^ ( / e^^^^*dU ) A{\fd\. 
i=i i=i 



Z„7r"("-i)/2 

Because of the Harish-Chandra, Itzykson-Zuber integral |151 116j 

det fe^J^* 



n— 1 



e^^^^*dU= I Hi! 

j=0 



A(a)A(A) ' 



we obtain that 



Pn{z) 



-!(n-=oJ 



Z„7r"("-i)/2 
We expand the determinant 



" " det (p^o^k^ 

n(.-A,)ne-'"^''^^^A(AMA. 

i=i i=i 



(2.5) 



det (^e"^^*) = ^ (-1)<^ J^e'^^^-O) 



where is the symmetric group. Hence 

,^ . 2 



-'(n-=oj! 



Z„^"("-i)/2A(a) 



E (-1)'' / - n ^"""^'^^ n e"^^"-'^' A(A)dA. 



We make a change of variables A^- = Xtj{i)- Then (— l)'^A(A) = A(A'), hence in the sum over S^., 
we have n! equal terms and, by dropping the prime, we obtain (|2.1() with constant 



Zn — 



7r"("-i)/2A(a) 



(2.6) 



This proves part (a). 
Observe that 



n(^-A,)A(A) =A(A,z) 



1 Ai 

1 \n 
1 Z 



A? 
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and 



-(y(Ai)-aiAi) _^^g_{V/(Ai)-aiAi) ... yr^-(V(Xi)-aiXi) 



g-(y(A„)-a„A„) ^^g-{V'{A„)-a„A„) 
1 Z 



\n„-(V{\n)-anXn) 



Then ()2.2() follows immediately from this and H2.1() . This proves part (b). 
From ((T^ it follows that 



mio mil • • • rriin 

mnO rUnl ■ ■ ■ rUnn 

rujo rriji ■ ■ ■ rrijn 



0, 



for every j = 1, . . . ,n. This proves (|2.3|) . To prove uniqueness of P„ satisfying (|2.3j) . observe 
that by equating the coefficients of in H2.2|) we obtain that 



mio mil • • • mi „_i 



mno m„i • • • m„,„_i 



(2.7) 



Let Pn{x) = + pn-ix"' ^ + ■ ■ ■ + Po and set p = (po ■ ■ ■ Pn-i) • Then the equations ()2.3|) 
are written in terms of the vector p as 



Mp =-m, M = ("ijfc)j=i,...,„.fc=o,...,n-i ' = irnjn)j=i,...,n- 

By ()2.7I) . det M / 0, hence p and therefore P„ is unique. 

Proposition 12.11 can be extended to the case of multiple ajs as follows. 



□ 



Proposition 2.2. Suppose A has distinct eigenvalues ai, i = 1, . . . ,p with respective multiplic- 
ities Hi so that ni + ■ ■ ■ + Up = n. Let n^*^ = ni + • • • + rij and n^^^ = 0. Define 

Wj{x) = 2;'>'.-ie-(^(^)-''^^), j = l,...,n, 

where i = ij is such that n^^~^^ < J ^ n^^^ and dj = j — n*-*"^^. Then the following hold. 

(a) There is a constant Z„ > such that 



" .7 = 1 .7 = 1 



(2. 



(b) Let 



fnjk = J X Wj{x)dx. 
Then we have the determinantal formula 

mio mil • ■ ■ fnir 



Pn{z) 



1 



m„o m„i 
1 z 



(2.9) 
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(c) Fori= l,...,p, 



Pn(2;)x^e-(^(^)-"'^)dx = 0, j = 0,...,ni-l, (2.10) 

and these equations uniquely determine the monic polynomial Pn- 
Proof. We write 

a = (ai,a2, • • ■ ,an) = ( ai, ■ ^ ,ai , 02, • • • ,ap-i,ap, . . . ,ap) 

ni times 



Oj's. In this limit we have that 



Apply formula (|2.5j) in the case when all aj = dj are different and take a limit to the multiple 

det (e'^J^fe) 



hm ■ 



det(Aj ^e"^^fe 



A(a) Ao(a)nr=ireri'fc! 



where dj is as in the statement of the proposition, and 



Thus, formula (|2.5() becomes 



Pn(2) 



det(Aj ^e°^^^ 



j=i j=i Ao(ajlli=illfc=i 



Then we continue as in the proof of Proposition 12.11 that is, we write 



f-A(A)dA. (2.11) 



det ( Xf'^e"^^^ 



o-eSn j=i 

and insert this into 1)2. 11[) to obtain a sum of n! equal terms, which leads to 1)2. 8() with 



U]=oJ 



(2.12) 



This proves part (a). 

Parts (b) and (c) follow from ()2.8() in the same way as parts (b) and (c) of Proposition 12.11 
followed from 1)2. 1|) . Note that in particular we have as in 1)2. 7() . 



mio mil • • • m-i „_i 



m„o m„i ■ • • m„^_i 



(2.13) 



□ 
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Remark: Formula (|2.8|) can be also written in the following form: 



where = (A„(i-i)_,_]^, . . . , A„(i)) and 

Zn = ZnUil . . .Upl. (2-15) 

When ^4 = 0, H2.14() reduces to the usual formula for Pn{z) with respect to the random matrix 
ensemble without external source. 

Corollary 2.3. Under the same assumptions as in Promsition \2.i\ we have that 

/oo 
P„(x)x"'e-(^(^)-"'^)(ix 7^ (2.16) 
-oo 

fori = l,... ,p. 

Proof. Let Pn+i be the average characteristic polynomial of an ensemble of (n + 1) x (n + 1) Her- 
mitian random matrices whose external source has the same eigenvalues as A plus an additional 
eigenvalue Oj. Then by part (c) of Proposition 12.21 we have that Pn+i is the unique monic poly- 
nomial that satisfies the relations ()2.3() with nj replaced by nj + l. If P„(x)x"'*e~'^^(^)~"'^^(ix 
would vanish, then Pn+i + Pn would satisfy these relations as well, which would contradict the 
uniqueness of Pn+i- □ 

Remark: The relations (|2.l'{j) can be viewed as multiple orthogonality conditions for the poly- 
nomial Pn- There are p weights e~^^^^^~"'^^\ j = l,...,p, and for each weight there are a 
number of orthogonality conditions, so that the total number of them is n. This point of view is 
especially useful in case A has only a small number of distinct eigenvalues. We will come back 
to this in Section 5 when we are considering the case of two distinct eigenvalues in detail. 

There is a considerable literature on multiple orthogonal polynomials (also called Hermite- 
Pade polynomials), see e.g. ^121^] and the references therein. 

3 Determinantal form of joint probability density function 

As in ProDosition l2.21 we assume that yl is a fixed Hermitian matrix whose eigenvalues oi, . . . , Op 
have respective multiplicities ni, . . . ,np, so that X]f=i ''^i — ^- S„ be the collection of 

functions 

$]„:={x^e"»" |i = l,...,p, j = 0,...,ni-l}. (3.1) 

We start with a lemma. 

Lemma 3.1. There exists a unique function Qn-i in the linear span of Tin such that 

/oo 
x^Qn^i{x)e''^^''Ux = forj = 0,...,n-2, (3.2) 
-oo 

and 

/oo 
x"-^Q„-i(x)e-^(^)dx = 1. (3.3) 
-oo 
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Proof. The conditions (|3.2|) and (|3.3|) give us n linear equations for the n coefficients of Qn-i 
with respect to the basis S„ with coefficient matrix 



mil rn2i 

V"^l,n-l "l2,n-l 



rrino \ 
rrini 



m. 



i.n— 1 / 



where nijk is as in part (b) of Proposition 12.21 This matrix is non-singular by (|2.13|) , so that 
the linear equations have a unique solution, and therefore Qn-i exists and is unique. □ 

For the rest of this section, we choose some ordering of the eigenvalues of A taking into 
account the multiplicities, say 

ai,a2, • • • (3.4) 

so that each Oj appears exactly rii times among the a's. For each k = 0,1,..., n, we can 
construct as in the previous section, but based on ai, . . . , a^. Thus is a monic polynomial 
of degree k such that 
roo 

e-(^(")-'^'")dx = 0, i = l,...,p, j = 0,...,ki-l, (3.5) 



/ 



where ki is the number of times that Oj appears among ai, . . . , a^. We also have tht P^ is the 
average characteristic polynomial of the ensemble of k x k Hermitian matrices with external 
source having eigenvalues Oj with multiplicity ki. 

For each A; = 1, . . . ,n we also have by Lemma l3.1l a function Qk~i from the linear span of 
the functions 

Sfc := {x^e'^'" \ i = l,...,p,j = 0,...,ki-l}. (3.6) 

such that 

roo 

x'Qfe_i(x)e-^(^)(ix = for i = 0,...,A;-2 (3.7) 



and 



-iQfe_i(x)e-^(^)dx = 1. 



(3.8) 



It follows from 1)3. 5|) . 1)3. 7|) . and (|3.8() that the P's and Q's are a biorthogonal system in the sense 
that 



Pj(x)Qfc(x)e-^(")dx = 5,fc, 



iov j, k = 0, . . . ,n — 1. 



(3.9) 



This property explains why we used Qk-i for the function that satisfies (|3.7|) and (|3.8j) (and not 

Qk). 

We now introduce the kernel Kn. 
Definition: With the polynomials Pk and the functions Qk introduced above, we define 

n-l 

Kn{x,y) = e'l(^(^)+^('^)) j;Pfc(x)Qfe(y). (3.10) 

k=0 

Note that the P's and the Q's depend on the specific ordering (|3.4jl that we choose for the 

eigenvalues of A. However, it will turn out that Kn does not depend on this ordering. 

Because of the biorthogonality property 1)3. 9|) it is easy to see from the definition ()3.10|) that 
we have 



Kn{x, x)dx = n 



(3.11) 
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and the reproducing kernel property 

Kn{x,y)Kn{y,z)dy = Kn{x,z). (3.12) 



The following is the main theorem of this paper. 
Theorem 3.2. The joint probability density function on eigenvalues has the determinantal form 

^det{Kn{Xj,Xk))i<j,k<n (3.13) 
The m-point correlation function has the form 

-Rm(Al, . . . , \m) = det(K„(Aj, \k))l<j,k<m (3-14) 

Proof. Any joint probability density function of the form (|3.13jl with a kernel Kn satisfying 
1)3. 11() and (|3.12|) leads to m-point correlation functions of the form ()3.14|1 . So it suffices to 
prove that (|3.13|) is the joint probability density function of the eigenvalues. 
For each j, we define 

Wj{x) = x'^^-^e"''' (3.15) 

if Oj = aj and Oj appears dj times in the sequence ai,. . . ,aj. Note that the functions (|3.15j) 
differ from the functions wj used in Proposition 12.21 in two respects. First there is an extra 
factor e"^*-^^ in Proposition 12. 2( and second we used a specific ordering of the eigenvalues of A 
in Proposition 12.21 (which only amounts to a renumbering). 

A similar calculation as that leading to (|2.5j) in the proof of Proposition 12.11 shows that the 
joint probability density of eigenvalues is proportional to 

n 

II e-^(^^) det{wi{Xj))i<i,,<nA{\). 
i=i 

Since Qi-i is a linear combination wi, . . . ,Wi we can take appropriate row combinations to 
find that 

det{wi{Xj))i<ij<n oc det((5j_i(Aj))i<jj<„. 
We write A (A) as a Vandermonde determinant which we similarly rewrite as 

A(A) = det(Pi_i(Afc))i<j,fc<„. 

Thus the joint probability density of eigenvalues is proportional to 

det (e-^^(^^)Q._i(A,)) det (e-^^(^'=)p,_i(Afc)) 

Taking the transpose of the matrix in the first determinant, and then using the multiplicative 
property of determinants, we find that the joint probability density is equal to 

cdet(ii:„(Aj, Afc)) 

l<j,k<n 

for some constant c, which should be such that the integral with respect to dXi ■ ■ ■ dXn is 1. 
Because of the properties (|3.11|) and (|3.12|) this is so for c = ^ and the theorem is proved. □ 

Remark: Renumbering the eigenvalues oi, 02, . . . On leads to the same kernel Kn but to different 
Pfc and Qk- 
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4 Special form of the kernel in case of two eigenvalues 

In this section we assume we have only two distinct eigenvalues ai and 02 with multiplicities ni 
and n2, respectively. We order the eigenvalues ai, 02, • • • , a„ in some arbitrary way {aj appear 
rij times in the sequence), but for convenience we assume that 

an-i = ai an = 02- (4.1) 

We also put 

a„+i = ai a„+2 = 02- (4.2) 

As in the preceding section, we have polynomials and functions Qk for every k = 0, . . . ,n — 1 
such that 

71—1 

K^{x,y) = e-"2(y(-)+y(y))J2Pk{x)Qk{y). 

k=0 

It is our aim in this section to simplify this expression. The formula we will find is an analogue 
of the well-known Christoffel-Darboux formula for orthogonal polynomials. 

To present the formula we are going to use multi-index notation. For non-negative integers 
ki and k2, we use Pk^M to denote the monic polynomial of degree ki + k2 having kj orthogonality 
relations with respect to the weight 

t^^.(x) = e-(^(^)-'^^-^\ j = l,2. 

Thus 



J —< 



Pki,k2ix)x^Wj{x)dx = 0, i = 0,... ,kj -1, j = 1,2. (4.3) 

-00 

The polynomial PkiM is called a multiple orthogonal polynomial of type II, see e.g. ^El- We 
also define 

Qk.Mi^) = Afc,,fc,(x)e'^i^' + Sfc,,fc,(x)e'^2^' (4.4) 
where the degree of Ak^^k2 is /ci — 1, the degree of -6^1,^2 is ^2 — 1 and 

f 0, j = 0,...,/ci + /c2-2, 
x^Qk,M{x)e-^^''^dx= I (4.5) 

[1 j = ki + k2-l. 

The polynomials Ak^^k2 and Bk^^k2 are called multiple orthgonal polynomials of type I, see 0I2]- 
For each pair (^1,^:2) of non-negative integers, the polynomials PkiM-- -^kiMi -^^1,^2 exist 
and are uniquely defined by their degree requirements and the relations (|4.3() . 1)4. 4() . and ()4.5() . 
We can express P^ and Qfc in this new notation as 

Pk = Pki,k2, and Qk-i = QkiM 

provided aj appears kj times among the numbers qi,... ,0^ (for j = 1,2). In particular, we 
have because of our assumptions 1)4. 1|) and 1)4. 2|) 

Pn — Pni,n2i Pn—1 — Pni,n2 — lj Pn—2 — Pni — l,n2 — l- (4-6) 

and 

Qn~l — Qn\,n2i Qn — Qn\+l,n2i Qn+1 — Qn\+l,n2+l' (4-7) 

We also need the numbers 

^klk2= Pk,M{^)x''^Wj{x)dx, j = l,2, (4.8) 



9 



which are non-zero, cf. H2.16|l . For later use we note that 

Pki,k2i^)i^ki+l,k2i^)^'i-i^) + Bki+i,k2ix)w2{x))dx 
Pfci.fca (2;)^fci+i,/c2 {x)wi {x)dx 



so that 



Similarly, 



(leading coefficient of Ak^+i^k2) x ^kiM 



leading coefficient of Afc^+i k2 = ~it\ — • (4-9) 



leading coefficient of Bk. fc^+i = —pr^. — • (4-10) 

h)' 



k\,k2 

It also follows from and dlHH) that h^f!^^^^ / for j = 1, 2. 
Then we can state the following theorem. 

Theorem 4.1. With the notation introduced above, we have 

I ^ni ,n2 



ri\ Pni~l,n2{x)Qni+l,n2 



1 

"'rai-l,n2 
^(2) 

~Pn\,n2 — l (•^)Qni ,n2+l (y) (4.11) 



"■ni,n2-l 



The proof of the theorem needs some preparation. We start working again with the P^'s and 
Qfs (single index) as before. For each j and /c, we put 



Co- 



jk = j xPk{x)Q,{x)e-^^''^dx. 
The coefficients Cjk appear in the expansion 



fc+i 

xPk{x) = Y,CjkPj{x), (4.12) 



j=0 

since by the biorthogonality relation we have indeed 



Cjk= I xPk{x)Qj{x)e ^'^''Ux. 
Similarly, we have for j = 0, . . . , n — 1, 

xQjix) = ^ CjkQkix). (4.13) 



n+l 



k=0 



Note that by adding the two values a^+i and an+2 as we did in ()4.2|) . we have this expansion 
for every j < n — 1. 
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Lemma 4.2. (a) If j > k + 2 then Cjk = 0. 

(b) If k > j + 3 and if both oi and 02 appear at least once among aj+2,aj+3, • • • then 
Cjk = 0. 

Proof, (a) We have that 

J P(x)Qj(x)e-^(^)(ix = 

for every polynomial P of degree < j — 1. Since xPk is a polynomial of degree k + 1, it follows 
that Cjk = Oiik + l<j — 1. This proves part (a). 

(b) Let k and j be such that the conditions of part (b) are satisfied. Suppose that ai 
appears ki times among ai, . . . , Uk, and ji times among ai, . . . , aj+i. We put k2 = k — ki 
and j2 = J ' + 1 — ji- It follows from the assumptions that ji < ki and J2 < /c2- Then 
Qj{x) = QjijaC^) = ^iij2(^)^"^^ + -^ii j2(^)^"^^ where ^jija degree ji — 1 and -Bjua 
has degree j2 — 1- It follows that 



and xAj^j^{x) has degree ji < fei — 1 and xBj-^^j^{x) has degree j2 < A;2 — 1. Thus, by the 
multiple orthogonality property of Pk = PkiM-: have 

Pfc(x)xQj(x)e-^(^)(ix = 0. 

This proves part (b). □ 

We also need the following relations between near-by P's and Q's. 
Lemma 4.3. We have 

p _ ''ni-l.na-l _p n 

-'ni — 1,712 — 1 (X) l-'ni— 1,712 -'711,712 — 1^ 

^ni-l,n2 
4?-ln2-l 

= ^ (-P?ii-l,n2 ~ -fni,n2-l) (4-14) 

^711,712-1 



and 



Qni+l,n2+l = ^Y) (Q?ii,n2+1 ~ Q?ii+l,7i2; 

^ni ,712+1 



(Q?ll, 712 + 1 Q7ll+l,712 ) (4.15) 



"'ni+l,n2 



Proof. Since P„^_i^„2 ^^^d P„^^„2-i ^^e both monic polynomials of degree n, their difference is 
a polynomial of degree < n — 1. Since this difference has Uj — 1 orthogonality conditions with 
respect to wj for j = 1, 2, it must be a multiple of Pni-i,n2-i- Thus 

-f 711 — 1,712 Pni,n2 — 1 — T-fni — 1,712 — 1 

for some 7. Integrating this equation with respect to and x^^~^W2{x), we get 

= 7/ii^i^_i,„2-i' and -h^^ln^-i = lh^S-i,n2-v respectively This gives (gH. 
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Next we note that we have 



and also 



J {Qni,n2+1 - Qni+l,n2)e ^''''^ dx = - = 0, j = 0, . . . , m + 712 - 1 

I x"i+"2(Q„,,„,+i - Qn,+i,„,)e-^(")dx = 1-1 = 0. 



Since Qni,n2+i{x) — Qni+i,n2{x) = A{x)e"''^^ + B {x)e"'^^ where A has degree ni and B has degree 
n2, it follows that Qni,n2+i - Qni+i,n2 IS a multiple of Qm+i.na+i, say 

Qni,n2 +1- Q ni+l,n2 — /3Qni+l, 712+1- 

This means for the j4-polynomials that 

^ni,n2+l ^ni+l,rt2 ~ P-^ni+l,n2+l 

and looking at the leading coefficient (= coefficient of x"^) we get 

^ _ leading coefficient of Ani+i,n2 _ ^i?,n2+i 
leading coefficient of /^W^^ 

where we used (|4.9|) . We also get by considering the B-polynomials that 

^ _ leading coefficient of Bn^^n2+i _ ^ni+i.n2 
leading coefficient of h^ln2 

because of ^JU^. This proves (|ITT^ . □ 

Now we are ready for the proof of Theorem 14.11 

Proof. We note that (x — y)e2^^^^^^'^^y^^ Kn{x,y) has a telescoping character. Indeed we have 
by Km and Km . 

n—l n—l n—l 

{x -y)^Pk{x)Qk{y) = ^xPk{x)Qk{y) -^yPj{x)Qj{y) 
fc=o fc=o i=o 

n— 1 fc+l n— ln+1 

= - '^'^CjkPj{x)Qk{y) 

k=0 j=0 j=0 fc=0 

n—l n—l 

= Cn,n~lPnix)Qn~l{y) - CjnPj {x)Qn{y) - Cj ^n+lPj{x)Qn+l{y) ■ 

j=0 j=0 

Now observe that Cn,n-i = !> and that Cjn = for j = 0, ... ,n — 3 and Cj_n+i = for j = 
0, . . . , n — 2, which follows from Lemma 14.21 Thus 

-Cn-2,n-Pn-2(a;)Qn(y) 
-Cn-l,n-Pn-l(2;)Qn(y) 

-Cn-l,n+l-fn-l(a^)Qn+l(y)- (4.16) 
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In formula (|4.16j) we have reduced the n-term expression to four terms, which is already quite 
nice. However, we want to reduce to three terms only. Changing back to multi-index notation 
and using (|4.6|1 and ()4.7|) . we see that ()4.16|1 leads to 

C-n—2,nPn\ — l,n2 — liS^)Qn\+l,n2 iv) 
C-n—l,nPn\ ,n2 — 1 (''^)Qni+l,n2 (?/) 

{x)Qni+l,n2+l{y) (4-17) 

Comparing (|lTrj) and we see that we need to get rid of P„j_i^„2_i and Qni+i,n2+i- 

This can be done using the following relations between near-by P's and Q's. 

Our next task is to express the recurrence coefficients Cn-2,n, Cn-i^n-, and c„_i^„+i that appear 
in (|4.17|) in terms of the /i-numbers. This is rather straightforward from the definition. Indeed, 
we have 



Cn-2,1 



X-P^i ,n2 (■^)Qni,n2 — 1 ('^)^ ^ ^ dx 

Pn-i,n2{x) {xAn^^n2-i{x)wi{x) + xBn^^n2-l{x)w2{x)) dx 
Pni,n2 {x)x-^ni,n2 — l {x)wi {x)dx 



(leading coefficient of An-^^n2-i) ^ ^nl 
h 



'ni,n2 

(1) 

^ni,n2 



"'ni-l,n2-l 



(4.18) 



where we used (|4.9|) . and similarly, 



"■?ii,n2 "-ni,n2 



Cn-l,n - + (2) ' ^^-^^^ 

"'ni-l,n2 "'ni,n2-l 



and 



_ "-ni+l,n2 / - r,„x 

Cn-l,n+l — — • l4.2Uj 

^7ii,n2-l 

Now we plug ah formulas (I^TTH) . (jlHSJ), (EUHl), <FT^ and into (|CTfl) . Straightforward 

calculations then lead to (14.111). □ 



5 Riemann-Hilbert problem 

We use the notation of Section 4. 

The Christoffel-Darboux formula (|4.11|) can be expressed in terms of the solution of a 
Riemann-Hilbert problem that was given by Van Assche, Geronimo, and Kuijlaars JHI to charac- 
terize the multiple orthogonal polynomials, and which generalizes the Riemann-Hilbert problem 
for orthogonal polynomials due to Fokas, Its, and Kitaev [T^. The Rieman-Hilbert problem is 
to find y : C \ M ^ C^''^ such that 

• y is analytic on C \ M, 
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for x E M, we have 



Y+{x) =y_(x) (01 
.0 1 



(5.1) 



where Y^{x) (Y^{x)) denotes the hmit of Y{z) as z — > x from the upper (lower) half-plane, 
as z — > oo, we have 



Y{z) = 7 + 





























(5.2) 



where / denotes the 3x3 identity matrix. 
In it was shown showed that there is a unique solution 

,n2Wl) C{Pn^^n2'W2) \ 



Y 



ClPni-l,n2 CiC{Pn^^l,n2Wl) CiC(P„^_i,„2 li;2 
\c2Pnun2-l C2C{Pn^^n2-l'Wl) C2C {Pn^^n2-l'W2) J 



(5.3) 



with constants 



,(1) 



ci = -27ri ^/i„,_i,„2j ' and C2 = 
and where Cf denotes the Cauchy transform of /, i.e., 

f{s) 



-2m K 



A2) 



(5.4) 



Cf{z) = / 



s — z 



■ds. 



The multiple orthogonal polynomials of type I Ani,n2J Bni,n2 have a Riemann-Hilbert char- 
acterization as well. We seek X : C \ M — > C^^'^ such that 



X is analytic on C \ ^ 
for X € M, we have 



as z ^ 00, we have 



1 0^ 

X+{x) = X_(x) I -'u;i(x) 1 

-W2{x) 1^ 



X{z) ={l + 



(5.5) 




— n 

























(5.6) 



The solution to this Riemann-Hilbert problem fT^ is 

/ -2'KiC{An^^n2Wl + BnT^^n2W2) '^'^i^m ,n2 27riB„^^„2 \ 



X 



-fclC(A„j+i^„2^1 + ^rji+l,n2l^2) ^l^ni+l,n2 ^l7?ni+l,n2 
\— A;2C(A„j_„2_|_iU'2 + -Bni,n2+1'"^2) ^2^ni,n2+l ^2-Bni,n2+l/ 



(5.7) 
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where ^ 

leading coefficient of ylrii+i,n2 ^"i'"-2' ^^'^^ 
and ^ 

^ leading coefficient of i?ni,n2+i "^'"^ ^ ^ 

It is easy to see that 

X = y~* (inverse transpose) (5.10) 

see also P^ . 

Now we form the product Y~^{y)Y{x) = X^{y)Y{x) and we compute the 21-entry using 

(E3I), (EH), (EU, (EH), 

(Pni,n2 (^) 
CiP„i_i,„2(x) 
C2-Pni,n2-l(a^)/ 

/ /lff,n2 /lSn2 ' 

= 27ri I P„j^„2(x)^ni ,712(2/) -fni -1,712 (^)^ni +1,712(2/) (2)~^ -fni ,n2-l (3^)^?ii ,n2+l (y) 

V ^ni-l,n2 ^ni,n2-l / 

We get a similar expression for the 31-entry [Y~^{y)Y{x)]3i, but with the 5-polynomials instead 
of the ^-polynomials. Then it follows that we can rewrite the Christoffel-Darboux formula (|4.1H) 
as 

K„(.,y) = ,-^(V».n»))^-'y-(!/)r(.)]..+e-».y-'(riV»]3i 

2-Ki[x — y) 

which is a compact form for the kernel in terms of the solution of the Riemann-Hilbert problem. 
We expect that the Riemann-Hilbert problem for Y is tractable to asymptotic analysis using 
the methods of or [TTllT^ rn]. 
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